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\ The field-theory model is proposed to study the electronic states near the Fermi energy in 

. spheroidal fullerenes. The low energy electronic wavefunctions obey a two-dimensional Dirac equa- 

tion on a spheroid with two kinds of gauge fiuxes taken into account. The first one is so-called K 
spin fiux which describes the exchange of two different Dirac spinors in the presence of a conical 
singularity. The second fiux (included in a form of the Dirac monopole field) is a variant of the ef- 
fective field approximation for elastic fiow due to twelve disclination defects through the surface of a 
spheroid. We consider the case of a slightly elliptically deformed sphere which allows us to apply the 
perturbation scheme. It is shown exactly how a small deformation of spherical fullerenes provokes 
an appearance of fine structure in the electronic energy spectrum as compared to the spherical case. 
^ ■ In particular, two quasi-zero modes in addition to the true zero mode are predicted to emerge in 

lyj ' spheroidal fullerenes. An additional 'hyperfine' splitting of the levels (except the quasi-zero-mode 

states) is found. 



O 



I 



X 



I. INTRODUCTION 



^ ■ Geometry, topological defects and the peculiarity of graphene lattice have a pronounced effect on the electronic 
O ' structure of fullerene molecules. The most extensively studied Cgo molecule is an example of a spherical fullerene 
nicknamed a "soccer ball" [1]. The family of icosahedral spherical fullerenes is described by the formula 2Q{n'^ + nl + P) 
with integer n and I. Other fullerenes are either slightly (as C70) or remarkably deformed and their general nickname 
. is a "rugby ball". The electronic structure of C70 cluster of D^^h geometry has been studied in [2] by using the 
^ local-density approximation in the density-functional theory. It was clearly shown that the spheroidal geometry is 
of a decisive importance for observed peculiarities in electronic states of the C70. What is important, the calculated 
energy levels were found to be in good agreement with photoemission experiments. 
\ Among different theoretical approaches to describe graphene compositions such as nanotubes, cones, nanohorns 
' and fullerenes the continuum models play a special role. Indeed, the continuum models allow us to describe some 
\^ ■ integral features of similar carbon structures. For instance, in the case of the spherical fullerenes the obtained 
I results can be extended on the whole 60n^ family. On the other hand, the continuum description gives an important 
information about long-distance physics which is difficult to get within other approaches like, e.g., widely- used tight- 
binding models. In addition, the continuum models are useful for clarification of pure topological effects such as an 
appearance of the Aharonov-Bohm phase and anomalous Landau levels (see, e.g., Refs. [3-5]). 

It should be noted that the formulation of any continuum model begins from the description of a graphite sheet 
(graphene) (see, e.g., [6-11]). The models which involve the inffuence of both geometry and topological defects on 
Q ] the electronic structure by boundary conditions were developed in Refs. [12,13]. A different variant of the continuum 
Q . model within the effective-mass description for fullerene near one pentagonal defect was suggested in Refs. [14,15]. 
' There are attempts to describe electronic structure of the graphene compositions by the Dirac- Weyl equations on 
the curved surfaces where the structure of the graphite with pentagonal rings taken into account is imposed with 
the help of the so-called K spin contribution [16,3]. In the case of spherical geometry, this is realized by introducing 
^ , an effective field due to magnetic monopole placed at the center of a sphere [16]. In the present paper we use a 
similar approach. The difference is that we consider also the elastic contribution. Indeed, pentagonal rings being 
the disclination defects are the sources of additional strains in the hexagonal lattice. Moreover, the disclinations 
are topological defects and the elastic flow due to a disclination is determined by the topological Frank index. For 
this reason, this contribution exists even within the so-called " inextensional" limit which is rather commonly used in 
description of graphene compositions (see, e.g., [17]). 

Recently, in the framework of a continuum approach the exact analytical solution for the low energy electronic 
states in icosahedral spherical fullerenes has been found [18]. The case of elliptically deformed fullerenes was studied 
in [19] where some numerical results were presented. In this paper, we suggest a similar to [18] model with the 
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Dirac monopole instead of 't Hooft-Polyakov monopole for describing elastic fields. We consider a slightly elliptically 
deformed sphere with the eccentricity of the spheroid e ^ 1. In this case, by analogy with [20] we use the spherical 
representation for the eigenstates, with the slight asphericity considered as a perturbation. This allows us to find 
explicitly the low- lying electronic spectrum for spheroidal fuUerenes. 



II. THE MODEL 



Our studies cover spherical molecules which are slightly elliptically deformed. Let us start with introducing 
spheroidal coordinates and writing down the Dirac operator for free massless fermions on the Riemannian spheroid 
S^. The Euler's theorem for graphene requires the presence of twelve pentagons to get the closed molecule. In a 
spirit of a continuum description we will extend the Dirac operator by introducing the Dirac monopole field inside 
the spheroid to simulate the elastic vortices due to twelve pentagonal defects. The K spin flux which describes the 
exchange of two different Dirac spinors in the presence of a conical singularity will also be included in a form of 
t'Hooft-Polyakov monopole. 

To incorporate fermions on the curved background, we need a set of orthonormal frames {ea}, which yield the same 
metric, g^v, related to each other by the local 50(2) rotation, 

ea^e'^= kiep, G SO{2). 

It then follows that 5^1/ = e^e^^a/s where is the zweibein, with the orthonormal frame indices being a, /? = {1,2}, 
and the coordinate indices n,^ = {1, 2}. As usual, to ensure that physical observables are independent of a particular 
choice of the zweibein fields, a local so(2) valued gauge field is to be introduced. The gauge field of the local 
Lorentz group is known as a spin connection. The Dirac equation on a surface S in the presence of the abelian 
magnetic monopole field is written as [21] 



where = + fl^ with 



n"e/[V^-z^^]V = i?^, (1) 



% = l^Vhc.,i0], (2) 

being the spin connection term in the spinor representation. 
The elliptically deformed sphere or a spheroid 

X V z 

— + ^ + - = 1, 3 

may be parameterized by two spherical angles = ^, = 6 that are related to the Cartesian coordinates x, y, z as 
follows 

X = a sin 9 cos (p; y = a sinO sin <p; z = c cos 6. (4) 
The nonzero components of the metric tensor for spheroid are 

= sin^ 6; g$$ = c? cos^ Q -\- (? sin^ ^, (5) 
where a,c>0, 0<^<7r, O<0<27r. Accordingly, orthonormal frame on spheroid is 

ei = — ^ 62 = , de , (6) 

a sm6i Va^cos^i+c^sin^ 

and dual frame reads 



= a sin 9 dtj); f? = \/ cos^ + sin^ Q d9 . (7) 
A general representation for zweibeins is found to be 



e 



% = a sin6'; e^g = 0; e\ = 0; e^e = V a'^ cos'^ 9 + sin^ 9 . (8) 
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Notice that e**^ is the inverse of e°j^. The Riemannian connection with respect to the orthonormal frame is written 
as [22,23] 

de^ = -w 2 A = -^^^=^^^^= d(j) A e^, (9) 
V cos^ + sin'^ 9 

de^ = -uj\ A = 0. (10) 

Here A denotes the exterior product and d is the exterior derivative. From Eqs.(9) and (10) we get the Riemannian 
connection in the form 



12 (f COS 12 

^<P2 = -^<pi = / 2 2n , 2 =^ ' '^^2 = '^^1 = ^^^^ 

V a cos^ + sin ^ 



We assume that the eccentricity of the spheroid e = -^/jl — (c/a)^! is small enough. In this case, one can write down 
c = a + Sa where 5 (| (5 |<C 1) is a small dimensionless parameter characterizing the spheroidal deformation from the 
sphere. So, we can follow the perturbation scheme using 6 as the perturbation parameter. 
Within the perturbation scheme, the spin connection coefficients are rewritten as 

Lol2 = 003 6(1 -6 sin^e). (12) 

The inclusion of the K-spin connection can be performed by considering two Dirac spinors as the two components 
of an SU{2) color doublet ij) = (V'tV'i)"^ (see Ref. [16] for details). In this case, the interaction with the color magnetic 
fields in a form of nonabelian magnetic flux ('t Hooft'-Polyakov monopole) is responsible for the exchange of the two 
Dirac spinors which in turn corresponds to the interchange of Fermi points. As was shown in Ref. [24], the Dirac 
equation for can be reduced to two decoupled equations for tp^ and ipi, which include an Abelian monopole of 
opposite charge, g = ±3/2. 

The elastic flow through a surface due to a disclination has a vortex-like structure [14] and can be described by 

Abelian gauge field W^. Similarly to Ref. [16], wc replace the fields of twelve disclinations by the effective field of the 
magnetic monopole of charge G located at the center of the spheroid. The circulation of this field is determined by the 
Frank index, which is the topological characteristic of a disclination. Trying to avoid an extension of the group we will 
consider the case of the Dirac monopole. The another possibility is to introduce the non- Abelian 't Hooft-Polyakov 
monopole (see Ref. [18]). As is known, the vector potential around the Dirac monopole have singularities. To 
escape introducing singularities in the coordinate system let us divide the spheroid (similarly as it was done for the 
sphere in [25]) into two regions, Bn and Rs, and define a vector potential in Rn and a vector potential in 
Rs- Notice that A^ includes both gauge fields. One has 

Rn: < 6* < I + A; < (/)< 277, (13) 
Rs: I - A < 61 < tt; < < 27r, (14) 

Rns : |-A<6i<| + A;0<(A<27r, (overlap) (15) 

where A is chosen from the interval < A < 7r/2. In spheroidal coordinates (5) the only nonzero components of A^ 
are found to be 

A^ g cos 9(1 + 6 sin^ 6i) + G(l - cos 6i) - sin^ 9 cos 9, (16) 

Af w 5cos6»(l + ^sin2 6i) -G(l + cos6») -(5Gsin2 6icos6i, (17) 

where the terms of an order of 6^ and higher are dropped. In the overlapping region the potentials A^ and A^ are 
connected by the gauge transformation 

A^ = A^^+iSNsd^S^l, (18) 
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where 



Sns = e^'^f 

is the phase factor. The wavefunctions in the overlap Rns are connected by 

ipN = SnsiPs, (19) 

where V'jv and ips are the spinors in Rn and Rs, respectively. Since Sns niust be single- valued [26], 2G takes integer 
values. Notice that the total flux due to Dirac monopole in (16) and (17) is equal to 47rG. What is important, there 

is no contribution from the terms with S. In our case, the total flux describes a sum of elastic fluxes due to twelve 
disclinations, so that the total flux (the modulus of the total Frank vector) is equal to 12 x 7r/3 = 47r. Therefore, one 
obtains G = 1. 

Let us consider the region R^- The Dirac operator for spheroidal fullerenes with monopole fields inside the spheroid 
takes the following form: 

:D = r»o + '5:Di, (20) 

where 

Vo = ij2-{de + ^) + i^^id^ - iG{l - cose) - ig cosO), 
a 2 asmt^ 

is the Dirac operator for spherical fullerenes and 



Vi = -171 



1 



a sin 9 



cos 9 sin 



7172 +iig — G) sin^ 9 cos 9 



sin^ ' 



- ij2 de, 



is the perturbation part of the spheroidal Dirac operator. This operator is hermitian on the spheroid. Actually, we 
consider the Dirac operator on the spheroid as a perturbation of the Dirac operator on the sphere (cf. [20]). For this 
reason, V should be transformed into the operator which is hermitian on the sphere. After straightforward calculations 
one obtains that the form of the operator Vq remains the same while Vi is written as 

©1 =-^sin6l(i-2mcos6'). (21) 

Let us turn back to Eq.(l) and formulate the eigenvalue problem for Dirac operator on the spheroid. Since a surface 
is a two dimensional space, the Dirac matrices can be chosen to be the Pauli matrices, 71 = — cr2,72 = —ci- We can 
restrict our consideration to only one of spinors, say to Vt- By using the substitution 

we obtain the Dirac equations for the spinor functions uj and vj in the form 

where m = g — G. To first order in S, the square of the spheroidal Dirac operator is written as = {Vq + dt), where 
r = {Voi^i + i^ii^o)- In an explicit form 

9^9 1 . ^.^ 1 1 TO— 7 cos ^ (7 — TO cos ^)^ 

sm^ 4 4sm^^ sm^ ^ sm'' 6* 

and 

a^f = 2f - jx (0-3 - 6m) - 4m(TO - 0-3)0;^ - 2m£73. (25) 

Let us write the equation V'^^p = E^tp by using the appropriate substitution x = cos9 in Eqs.(24) and (25). One 
obtains 
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(26) 



where A = aE. When 6 = 0, one has the case of a sphere with magnetic monopole inside. In this case, the exact 
solution is known (see, e.g., [18]) 



_ ^C„(l-x)"(l + a;)/3p2".2/3(^) 



with the energy spectrum 



Here 



iX%r = {n+\j\ + l/2f-m\ 



(27) 



1 

"=2 



1 — m — - 
2 



m+- 



, = — 
2 



J +m- - 



(28) 



P^"''^^{x) and P^i^''^'^{x) are Jacobi polynomials of n-th order, and C„ and Cy are the normalization factors. These 
states are degenerate. For example, the degeneracy of the zero mode is equal to six. Therefore, we have to use the 
perturbation scheme for the degenerate energy levels [27]. As a result, the energy spectrum for spheroidal fullerenes 
is found to be 



(Af J2 = (n + lil + 1/2)2 _ ^2 ^ ^^2f + jL' + jL-' + _ 2m£>,-„). 



(29) 



Similarly to [18], the possible values oi j obey the condition jjj > ||m| + 1/2] for non-zero modes. Direct calculations 
show that 



j(6m2 + l/2) jm2(4m2 + 5) 



p{p + 1) 



p2(p+ 1)2 



(30) 



with p = n + (} + a, 



L'-' = 4m(m- l)|C„|2F„(2a,2/3)/„(2a,2/?) 
+4m(m + l)\CyfFn{2fi, 2i/)/„(2^, 2iy) 



(2a, 2/3) = 



(n + l){n + 2a+ l){n + 2(3+ l)(n + 2a + 2/3 + 1) 
(2p+l)(p+l)2(2p-n + 3) 

n(n + 2a) (n + 2/3) (n + 2a + 2/3) 



+ - 



(2p-l)p2(2p+l) 



(31) 
(32) 



/„(2a,2/3) = 



22a+2/3+ip(-^ + 2a + l)r(n + 2/3+1) 
n!(2p+ l)r(n + 2a + 2/3+ 1) ' 



(33) 



Djn 



T{n + 2„+ l)r(n + 2^ + 1) - (j^^) ' ^t;;:;::::;;^ r(n + 2a + l)r(n + 2/3 + 1) 
r(n + 2m + l)r(n + 2Z/ + 1) + (^^^) ' ^(n + 2a + l)r(n + 2/3+1) 



(34) 



for j > and 



5 



_ r(n + 2m + l)r(n + 2^ + 1) - (^jj^) ^tljl;::::^;^ r(n + 2a + l)r(n + 2/3 + 1) ^^^^ 

r(n + 2m + l)r(n + 2Z. + 1) + ( 1^1+^^+" ) ' r(n + 2a + l)r(n + 2/3+1) 

for j < 0. Finally, in the linear in 6 approximation, the low energy electronic spectrum of spheroidal fuUerenes takes 
the form 



with 



E^, = ^ViMTnmTn, E>, = ^^. = '^"-^ ^ (37) 

jn jn 

where ^ = fj, {v) and r] = P {a) for j > {j < 0), respectively. Here we came back to the energy variable E = X/a (in 
units of TiVp/a where Vp is the Fermi velocity). 

Notice that the non-diagonal matrix element of perturbation (e^'^|(V-'jn|r|V-'-jri)k "''^) turns out to be zero. There- 
fore, the states with opposite j do not mix and j remains a good quantum number as would be expected. As is 
seen from Eq. (30), for 5 = both eigenstates and eingenvalues are the same as for a sphere (cf. Ref. [18]). At the 
same time, Eqs. (30)- (35) show that the spheroidal deformation gives rise to an appearance of fine structure in the 
energy spectrum. The difference in energy between sublevels is found to be linear in 5 which resembles the Zeeman 
effect where the splitting energy is linear in magnetic field. In addition we found a further splitting of the states with 
opposite j. The splitting is weakly pronounced ('hyperfine' splitting) and entirely dictated by the topological defects. 
In other words, this splitting is pure topological in its origin. The magnitude of hyperfine splitting is determined by 
the topological " charge" m. In particular, for m = (no defects) E'-^ = and no splitting occurs. 

Table 1 shows all three contributions (in compliance with Eq. (36)) to the first energy level for various fullerenes 
with different morphologies. As is clearly seen, the first (double degenerate) level becomes shifted due to spheroidal 
deformation. Schematically, the structure of the first energy level is shown in Fig.l. The more delicate is the structure 
of the second level which is presented in Table 2. In this case, the initial (for 5 = 0) degeneracy of E^^ is equal to 
six. The spheroidal deformation provokes an appearance of three shifted double degenerate levels (fine structure) 
which, in turn, are splitted due to the presence of topological defects (see Fig. 2). Notice that the magnitude of both 
shifts and the splitting energy increases with j and m. It should be stressed that the experimental observation of 
such splitting would indicate the presence of topological defects while a measured magnitude of the splitting gives the 
value of the monopole charge. Notice that in our model we have two possible charges: 771 = 5— G = ±3/2— 1. Foi' 
comparison, the model suggested in [24,16] predicts m = ±3/2. Thus, the splitting energy could clarify the proper 
monopole charge due to topological defects. Notice that a similar small splittings of the spheroidal subshells due to 
D^h symmetry was discovered in Ref. [2]. 

Let us discuss briefly the zero-mode state. As stated above, for a spherical fullerene there is a zero-mode state with 
a sixfold degeneracy. For spheroidal fullerene, the zero-mode state also exists, however, its degeneracy is twofold. In 
addition, there appear two slightly shifted "quasi-zero" modes for j = ±1, ±2, each of them is also twofold degenerate 
(there is no hyperfine splitting for quasi-zero- mode states). Indeed, our study shows that in this case Djn is equal to 
one, so that Ajn in (30) becomes odd function of j. It should be stressed that the similar conclusions follow from our 
analysis based on spheroidal harmonics suggested in Ref. [28]. Thus, we confirm the finding in [2] that there are only 
up to twofold degenerate states in spheroidal fullerene. Generally, the spheroidal deformation gives rise to a reduction 
in degeneracy of energy levels in comparison with the spherical case. 



III. CONCLUSION 



We have considered the electronic states of spheroidal fullerenes provided the spheroidal deformation from the 
sphere is small enough. In this ease, the spherical representation is used for describing the eigenstates of the Dirac 
equation, with the slight asplierieity considered as a perturbation. The using of the perturbation scheme allows us to 
find the exact analytical solution of the problem. In particular, the energy spectrum of spheroidal fullerenes is found 
to possess the fine structure in comparison with the case of the spherical fullerenes. The energy between sublevels is 
found to be linear in the small distortion parameter 5 and is positive/negative for prolate/oblate spheroid, respectively. 
Another principal difference is the obtained hyperfine splitting of the energy levels. We have shown that this splitting 
is weakly pronounced and entirely dictated by the topological defects, that is it has a topological origin. Our finding 
confirms the results of [2] that there can exist only up to twofold degenerate states in the C70. 
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Notice that for spherical fullerenes {S = 0) our results agree with those found in [18]. It is interesting that the 
predictions of the continuum model for spherical fullerenes are in qualitative agreement with tight-binding calcula- 
tions [29-32]. In particular, the energy gap between the highest-occupied and lowest- unoccupied energy levels becomes 
more narrow as the size of fullerenes becomes larger. It is important to keep in mind, however, that the continuum 
model itself is correct for the description of the low-lying electronic states. In addition, the validity of the effective 
field approximation for the description of big fullerenes is also not clear yet. Actually, this approximation allows 
us to take into account the isotropic part of long-range defect fields. For bigger fullerenes, one has to consider the 
anisotropic part of the long-range fields, the influence of the short-range fields due to single disclinations as well as 
the multiple-shell structure. Therefore, we admit that the proper values of the energy levels can deviate from our 
estimations. 

Finally, the main results obtained for spheroidal fullerenes are the discovery of (i) fine structure with a specific shift 

of the electronic levels upwards, (ii) hyperfine splitting due to topological defects, and (iii) three twofold degenerate 
modes near the Fermi level with one of them being the true zero mode. In our opinion, these predictions are quite 
general for the fullerene family and are of interest for experimental studies. 

The work was supported in part by VEGA grant 2/6193/26 of the Slovak Academy of Sciences, by the Science 
and Technology Assistance Agency under contract No. APVT-51-027904 and by the Russian Foundation for Basic 
Research under Grant No. 05-02-17721. 
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TABLE I. The structure of the first energy level for various fuUcrcncs with different morphologies given in Ref. [33]: S — spherical, I — 
icosahedron faceted, TI - truncated icosahedron; except Cigo - the structure given in Ref. [34] and YO - given in Ref. [35]. b is average 
bond length, R (R = a) is average radius, SD is standard deviation from a perfect sphere (see Refs. [33,34]), so that S = SD/R. The 
hopping integral is taJien to be t = 2.5 eV and Vp = 3tb/2h . 



YO-C240 


b = I.45A R = 7.O3A SD = O.I7A 5 = 0.024 


j 




\EU\{rneV) 


\E^„\{meV) 


n = 1, m = 


1/2 


1 


1.89 


6.5 


2.5 






-1 


1.89 


6.5 


2.3 


n = 0, m = 


1/2 


2 


1.89 


26 


2.2 






-2 


1.89 


26 


1.8 


n = 0, m = 


-5/2 


3 


1.89 


4.6 


13 






-3 


1.89 


4.6 


18 



TABLE II. The structure of the second energy level for YO-C240. 



9 



FIG. 1. The schematic picture of the first positive electronic level i5*„ for spheroidal fuUerenes. 
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FIG. 2. The schematic picture of the second positive electronic level for spheroidal fullerenes. 
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